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ABSTRACT 

The present knowledge on hard processes in the context of heavy ion colhsions is 
reviewed, with particular emphasis on J/ip production. The p-A data on charmonia 
production from Fermilab experiments is shown to be in excelent agreement with 
the p-A data collected at CERN. The simultaneous analysis of all existing p-A data 
reaches a precision which allows us to rule out some preconceived ideas, setting a good 
frame against which the data collected with ion beams at CERN can be compared. 



Summary talk presented at the 12th International Conference on Ultra-Relativistic 
Nucleus-Nucleus Collisions (Quark Matter '96), Heidelberg, Germany, May 1996. 



1 



Hard probes in nucleus-nucleus collisions 

Carlos Lourengo 

CERN PPE, CH-1211 Geneva 23, Switzerland 

The present knowledge on hard processes in the context of heavy ion collisions is re- 
viewed, with particular emphasis on J/?/' production. The p-A data on charmonia produc- 
tion from Fermilab experiments is shown to be in excelent agreement with the p-A data 
collected at CERN. The simultaneous analysis of all existing p-A data reaches a precision 
which allows us to rule out some preconceived ideas, setting a good frame against which 
the data collected with ion beams at CERN can be compared. 

1. Introduction 

The first data on i/ifj production in nucleus-nucleus collisions were collected 10 years 
ago, at the CERN SPS, by the NA38 experiment, with the 200 GeV per nucleon oxygen 
beam. 

Although the experiment had been planned to search for thermal dimuons, a paper 
by Matsui and Satz redirected attention by concluding that "there appears to be no 
mechanism for ^/if) suppression in nuclear collisions except the formation of a deconfin- 
ing plasma, and if such a plasma is produced, there seems to be no way to avoid ^/ip 
suppression." 

In fact, already at QM '87, the NA38 collaboration presented evidence of a substantial 
decrease of the ratio T/'/Continuum, both from p-U to 0-U collisions and versus i?T within 
0-U [0. These observations were soon confirmed with S-U data 

Although the exciting interpretation of this "J/?/' suppression" as a signal of QGP 
formation was always present, some alternative explanations of a more "classical" nature 
were soon proposed. 

Several developments helped to clarify the interpretation of the data: 

• Concerning the absolute cross sections, Fermilab experiment E772 observed that 
^/tp production is, in p-A collisions, proportional to A°, with an a value smaller 
than the one previously reported by NA3 (0.94±0.03 [|]) and used in early NA38 
papers [^f. 

• Several measurements of B^f^a'^ in p-A collisions were obtained by NA38, leading 
to a much more complete and accurate reference in the study of the ion data. 

• Higher statistics in the S-U data set allowed to use the high mass Drell-Yan as 
the reference in the dependent (relative) studies, replacing ^/i/Cont by ip/DY. 
Furthermore, the isospin correction (~ 10 %) required when comparing the proton 
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to the (isoscalar) sulphur induced colhsions lead to a smaller "suppression" relative 
to what had been reported before, using ip/Coni. 

At the same time, information on the ipi became available, revealing a "suppression" 
much beyond the nuclear absorption expected from the p-A data. 

• The absolute cross sections of direct i/ip and ■?/'/ production were found to be far 
too high relative to the values calculated with the colour singlet model. Although 
this was known to be the case already from fixed target data the new (high px) 
data from CDF 0, clearly identifying the fraction of ^/ip coming from beauty and 
Xc decays, triggered a major change in the theoretical understanding of charmonia 
production . Contrary to the paradigm of the colour singlet model, a substantial 
fraction of the ^/ip and ipi cross section is now believed to come from binding of cc 
pairs produced in a colour octet state. 

Taking into account all these developments, some of them only available after QM '95, 
Kharzeev and Satz were able to explain the phenomenological fit made earlier by 
Gerschel and Hiifner |T0[, by associating the absorption cross section of around 6-7 mb 
to the ccg state crossing the nuclear matter, rather than to a fully formed J/'?/', for which 
they had previously calculated a much smaller cross section |TT| . 



The further suppression of the ipi observed in S-U collisions is explained [^] by the 
formation of a (confined) hadronic medium in the later stages of the collision. These 
"comovers" can easily break the loosely bound, and by then fully formed, ipi while leaving 
the unchanged, since confined gluons are too soft to overcome the ^/ip^s (~ 640 MeV) 
binding energy. 

This was the status of the charmonia front just before the results from the Pb data 



collected in December 1995 became available. These new results certainly open 
another very exciting chapter in the ^/tp saga, and the people attending QM '96 will 
remember the lively discussions held inside and outside of the conference room. 

Some of the proposed "explanations" are not accurate enough to be seriously considered. 
For instance, it is not reasonable to reproduce the Pb-Pb data at the price of missing the 
p-A and S-U points. Also, we cannot mix the ip/Coni values of S-U (presented for the 
last time at QM '91) with the ip/TiY values of Pb-Pb, the latter variable certainly being 



much better (remember the enhancement observed in the intermediate mass region [13 
which might affect the "continuum"). 

The precision and diversity of the measurements now available certainly deserves equally 
accurate theoretical descriptions. The field of ^/ip production in heavy ion collisions has 
reached maturity and QM '96 should be considered as a transition point: the time has 
come to build realistic theoretical models that reproduce the data with an accuracy at 
the few percent level. 

In this sense, it would be very helpful to develop a microscopic model (like RQMD, 
VENUS or HIJING) that would generate hard processes in a nuclear environment, in- 
cluding nuclear geometry effects, modifications of the parton distribution functions, etc. 
Of course, it would have to work also at SPS energies, not only for RHIC and LHC. 

In the remaining part of this paper, the data presently available is presented in detail, 
both the directly measured values and, where they are required, those obtained after 
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Figure 1. p-A Drell-Yan double differential cross section, per target nucleon. See text for 
details. 



correcting for ^/s dependences and phase space windows. Special emphasis is given to the 
p-A data, collected at FNAL and at CERN, enough by it self to rule out some ideas. 



2. Hard processes in p-A collisions 

2.1. Drell— Yan production 

Figure [l| presents the double differential Drell-Yan cross section, d'^a/drndxp, in 
nbxGeV^ /nucleon, versus the scaling variable a/t = m/y/s (using mass bins of 0.5 GeV), 
for different xp bins (0.05 wide). For presentation purposes, the points above (below) 
the xp = 0.375 data set have been successively multiplied (divided) by factors of 10. 
This figure merges Drell-Yan data from three p-A dimuon experiments: E772, E605 and 
NA3 lUHll]. The dat available as tables from the hepdata data base |17i. The 



beam energy was 800 GeV in the Fermilab experiments and 400 GeV in the case of NA3. 
Different targets were used: hydrogen (E772), copper (E605) and platinum (NA3). The 
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agreement between the three experiments is excelent. The curves on the figure are the 
result of a leading order calculation using the CTEQ2M set of (NLL) quark distribution 
functions [jl8[, multiplied by an appropriate K factor (around 1.5). The K factor is 
considerably reduced when the calculations are performed in NLO and becomes ~ 1 at 
NNLO 

The infiuence of the nuclear target, usually referred to as "the A-dependence" , is quite 
simple in the case of Drell-Yan production. The cross section just scales linearly with 
the number of target nucleons, = a^^ x A. A closer look reveals that we have to 
count the number of protons and neutrons separately, since the Drell-Yan cross section 
depends on the charge of the annihilating quarks. In fact, the precision reached in Drell- 
Yan measurements is such that it has revealed differences between the structure functions 
of the u and d sea in the nucleons pO|, a probable explanation for the violation of the 



Gottfried sum rule recently observed in deep-inelastic lepton scattering [BTI]. 




0.0 0.3 0.6 



1.0 2.0 3.0 
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Figure 2. Drell-Yan A-dependence expressed as a versus xt, m, x-p and px- Notice the 
fine scale used in the vertical axis. 



The dependence of Drell-Yan on A, for p-A collisions, can be expressed in terms of the 
a parameter, as done in figure H, or by the ratios between the dimuon yields per nucleon 
in heavy targets relative to a light target, say deuterium, as presented in table |l|. 
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Table 1 

Ratios of DY, ^/ip and ipi cross sections between p-A and p-^H, from E772 data. 



Target (A) 


DY 






C (12) 


1.003±0.0085 


0.851±0.0125 


0.855±0.0292 


Ca (40) 


0.995±0.0061 


0.806±0.0087 


0.750±0.0271 


Fe (56) 


0.990±0.0065 


0.756±0.0102 


0.722±0.0357 


W (184) 


0.986±0.0083 


0.619±0.0133 


0.623±0.0331 



2.2. i/ip production 

Table |1| also includes the values measured by the E772 experiment I^Sj for the J/^ and 
ipl resonances, first shown at the QM '90 conference. 

Before comparing the results of the CERN experiments, NA38 and NA51, with the 
E772 data, in what concerns p-A collisions, we should notice that the a of i/ip depends 
on x-pi as is shown in figure H W^- The E772 data on is also presented in table El 
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Figure 3. xp dependence of a and . 



Since the i/ip cross section becomes negligible at high xf, the average value of a'^, for 
positive xf, is around 0.91. In fact, the value = 0.920 ± 0.008 was published by E772 
in Ref. ^ but the E789 collaboration uses the value = 0.90 ± 0.02 in Ref. E3|. 

We move now to the p-A measurements done at the CERN/SPS by the NA38 and NA51 
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Table 2 

Values of versus xp, from E772 data. 



x-p 




0.15 


0.919±0.005 


0.25 


0.918±0.002 


0.35 


0.910±0.004 


0.45 


0.883±0.008 


0.55 


0.866±0.009 


0.65 


0.815±0.034 



collaborations. These data are certainly the best to study the (target) nuclear effects 
on J/^/; production as a reference for the study of the heavy ion data, since they were 
collected with the same dimuon spectrometer, therefore being in the same kinematical 
window. However, we should start by checking that there is compatibility between the 
p-A results of NA38/51 and the observations done at Fermilab. 

The p-A results of NA38 and NA51, in what concerns J/^/^ cross sections, are presented 
in table |[ These values correspond to the phase space domain (A) covered by the detector, 
< ?/* < 1 and I cos6cs\ < 0.5, y* and 6cs being the center of mass rapidity and the 
dimuon decay angle in the Collins-Sopper frame. 



Table 3 

Cross sections for ^/ip production in p-A collisions, times b.r. into muons, as measured 
by the NA38 and NA51 collaborations. See the text for details. 





Pbcam 




Bal/A 




(GeV) 


(nb/nucleus) 


(nb / nucleon) 


PP 


450 


4.41±0.52 


4.41±0.52 


P-2H 


450 


9.16±1.05 


4.58±0.52 


p-C 


450 


46.7±5.3 


3.89±0.44 


p-Al 


450 


93.8±10.7 


3.47±0.40 


p-Cu 


450 


216±24 


3.43±0.38 


p-Cu 


200 


104±23 


1.65±0.37 


p-W 


450 


566±65 


3.08±0.36 


p-W 


200 


263±18 


1.43±0.10 


p-U 


200 


326±69 


1.37±0.29 



In order to fit the ^/ip cross sections measured in p-A collisions by the NA38 and NA51 
collaborations to the A" form, we must first correct for the different beam energies used. 
Indeed, the data was collected either with 200 GeV or 450 GeV proton beams. Since the 
heavy ion data was collected at 200 GeV/nucleon (158 GeV in the case of the Pb beam) 
we take i/i = 19.4 GeV as the reference and rescale the other values. 

The a/s dependence of the ^/ip cross section (times b.r. into muons) is displayed in 
figure ^ basically taken from fig. 11 of Ref. ^ but with a few more points. To make this 
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figure, the "per nucleon" cross sections were evaluated using an a of 0.91 (it would look 
essentially the same with 0.92). 




Figure 4. Energy dependence of the J/^/^ cross section. 



Notice that we have included in this figure the values measured by NA51 and NA38 
(after converting them to the positive domain, see below). They are in excelent 
agreement among them selfs and fit very well within the other measurements. 

The solid line on the figure (the dotted line will be explained later) is 

^(^"^ f ^^ A 3.097V' • t,/ 1 /-.X 
\y s) = 37 X II — I , m nb/nucieon. (Ij 



We have used this function to rescale the 450 GeV values, the result being given in 
table ^. In the last column of the same table are the corresponding values in the xp > 
phase space domain. We have multiplied the measured values by a factor of 2 to account 
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Table 4 

Cross sections for J/^ production in p-A collisions, times b.r. into muons, after corrections. 
See the text for details. 





Pheum 


-Ba^/A (rescaled) 


B(Tt^>o/A (rescaled) 




(GeV) 


(nb / nucleon) 


(nb / nucleon) 


PP 


450 


2.11±0.25 


4.52±0.54 


p-^H 


450 


2.19±0.25 


4.69±0.54 


p-C 


450 


1.86±0.21 


3.98±0.45 


p-Al 


450 


1.66±0.19 


3.55±0.41 


p-Cu 


450 


1.64±0.18 


3.51±0.39 


p-Cu 


200 


1.65±0.37 


3.53±0.79 


p-W 


450 


1.47±0.17 


3.15±0.36 


p-W 


200 


1.43±0.10 


3.06±0.21 


p-U 


200 


1.37±0.29 


2.93±0.62 



for the cos6cs coverage (assuming a uniform decay angle distribution) and by a factor 
1.07 to convert the y* range into the positive xp window. 

From the rescaled values presented in table we can fit an A" function and extract 
a^. The result is 0.908±0.029, excluding the pp point (otherwise it is 0.919±0.021). 

This is in excelent agreement with the results from Fermilab E772/E789 experiments. 

2.3. ijjf production 

ipf production in p-A collisions is easier to study through the ratio 'ip'/ip, since the 
J/ip has well understood energy and A dependencies. The values of the ratio ratio 
of production cross-sections times branching ratios into muon pairs, measured in p-A 
collisions |^,^], are collected in table ^. 

Table 5 



Values of the ip'/ip ratio obtained in p-A collisions. 







Pbeam (GcV) 


Ba"^' 1 Ba^ (%) 


Exp.t 


p 


H2 


450 


1.69 ±0.03 


NA51 


p- 




450 


1.80 ±0.03 


NA51 


p- 


-C 


450 


1.90 ±0.13 


NA38 


p- 


-Al 


450 


1.36 ±0.35 


NA38 


p- 


-Cu 


450 


1.68 ±0.11 


NA38 


p 


-W 


450 


1.59 ±0.13 


NA38 


p- 


-W 


200 


1.80 ±0.17 


NA38 


p- 


-u 


200 


1.77 ±0.22 


NA38 


p- 


-p 


v/i = 63 


1.9 ±0.6 


ISR 


p 


-Li 


300 


1.88 ±0.26 ±0.05 


E705 


p 


Be 


400 


1.7 ±0.5 


E288 


p- 


-Si 


800 


1.65 ±0.20 


E771 


p- 


-Au 


800 


1.80 ±0.1 ±0.2 


E789 
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These values show that the ratio between tpl and ^/ip production cross sections is 
independent of the target nuclei and (at high energies) of i/i. In fact, we can even add 
values obtained with other beam particles ||2^,|26| (tt"*", tt", p, 7) without changing the 
picture. 

This observation is a clear indication that these charmonia bound states are actually 
only formed as such outside the target nucleus, at least when they are fast enough (at 
positive xp)- It would certainly be very interesting to see what happens when the cc state 
is slow enough to become a J/V' or tpl within nuclear matter Present (p-A) dimuon 
experiments are unable to detect the correspondingly slow decay muons but this might 
soon be overcome by an "inverse kinematics experiment", sending the SPS Pb beam on 
a light target. 




Figure 5. Target dependence of the ip' /%p ratio, in p-A collisions. 



The data of table ^ have been plotted versus A in figure ^ and versus y/s in figure |^. In 
these figures, some of the points were slightly displaced in the horizontal axis (otherwise 
a few points would be indistinguishable). 

The dotted curves in these plots are the result from linear fits to the data, the slopes 
being compatible with zero for both the A and the a/s cases, (—1.6 ± 4.2) x 10~^ and 
(-8.1 ± 79) X 10-^ respectively. 

Even if would be smaller than by only 0.01, that would imply a decrease of V'VV' 
from pp to p-U of 5.6 %, already quite difficult to accomodate by these measurements. 
Therefore, models which predict a substantially different A dependence for the J/^ and 
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Figure 6. Energy dependence of the ip'/tp ratio, in p-A collisions. 



ip/ resonances, in p-A collisions, are in clear contradiction with the data. 

Figure ^ also contains a second (dashed) curve, which represents the evolution with ^/s 
of the ip' /ip ratio expected if we assume the ipf cross section to depend on ^/s according 
to equation |T], with the parameter 3.097 GeV replaced by 3.685 GeV. The fact that this 
curve (normalised in the plot at ^/s = 19.4 GeV) completely fails to go through the 
data points, shows very well how incorrect such an assumption is, and reveals that this 
parameter should not be directly associated to the mass of the particle under study. In 
fact, the best fit to the points in figure § is provided by 31.5 x (1 - 4.4/v^)'-' (dotted 
line in that figure). 



3. Hard processes in nucleus- nucleus collisions 

Drell-Yan differential cross sections (corrected for detector specific acceptance and 
smearing effects) are not yet available for nucleus-nucleus collisions. Work in this di- 
rection is still in progress within the NA38/50 collaborations. Therefore, it is not yet 
possible to compare the data with NLO calculations (not available in the form of event 
generators). Although the nuclear effects on the parton distribution functions are ex- 
pected to be small, in the kinematical window of NA38/50, they are certainly non-zero 
and it would be very important to actually measure them rather than just neglect their 
existence. For the moment, the fact that the same K factor (relative to a certain LO 
calculation) is required by the p-A and ion data is the best argument to say that 



Drell-Yan is a good reference in the studies of charmonia production by ion collisions. 
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The J/^ cross sections measured in nucleus-nucleus collisions, by the NA38 and NA50 
collaborations, are collected in table |, before and after ^/s and phase space corrections. 



Table 6 

J/i/j cross sections, times b.r. into muons, measured by the NA38 and NA50 collaborations 



in B-A collisions. 





Pbeam 


Bal 


Bal/{BxA) (rescaled) 


i?cr^p>o/(BxA) (rescaled) 




(GeV) 


(/ib) 


(nb/nucleon) 


(nb / nucleon) 


0-Cu 


200 


1.26±0.13 


1.25±0.13 


2.68±0.28 


0-U 


200 


4.42±0.46 


1.16±0.12 


2.48±0.26 


S-U 


200 


7.69±0.76 


l.OliO.lO 


2.16±0.21 


Pb-Pb 


158 


19.0±1.4 


0.59±0.04 


1.26±0.09 



The light ion collisions are in excelent agreement with the extrapolation of the proton 
data "nuclear absorption" curve. Indeed, if we include the oxygen and sulphur points in 
the fit to a (BxA)" function, the value of a stays the same (0.911±0.016) as before. This 
smooth transition from p-A to O, S-A collisions is illustrated in figure |^. 
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Figure 7. J/ip production cross sections (times b.r.) versus AxB. The straight line, in 
this log-log plot, corresponds to a = 0.908, the dotted part being an extrapolation from 
the fitted region. 
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Contrary to what happens in the hght ion region, the Pb-Pb point is seen to be away 
from the "expected" production cross section, which is 50 % higher than the measured 
value. This is a strong indication that something else besides "nuclear absorption" is 
happening to the ^/if) in Pb-Pb collisions. 

There are more precise ways of checking whether nuclear absorption alone can explain 
the observed evolution of J/?/' production yields. On the vertical axis we can use the 
ratio of J/?/' to Drell-Yan yields to reduce the systematic errors resulting from absolute 
normalisations. To explore the dependence of ^/if) production on the centrality of the ion 
collisions, we can subdivide the event samples according to some global variable, like E't 
or EzBC, sensitive to the geometry of the collision. In fact, the best variable within the 
"nuclear absorption" framework is the amount of nuclear matter crossed by the pre- J/?/' 
state on its way out, the famous L introduced by Gerschel and Hiifner [llO[ . 

Figure |^ attempts to merge in a single plot the variation with L of both the J/ip and 




o" 

A 

LL 

2^10 - 



41 



© p-p, d (450 GeV) (NA51 

ft p-A (450 GeV) (NA38) (* 

★ p-A (200 GeV) (NA38) 

A 0-Cu, U (200 A GeV) (NA38) 

□ 



(*) rescaled to 200 GeV 



S-U (200 A GeV) (NA38) ' 1 

Pb-Pb (158 A GeV) (NA50) (*) 



2 4 6 8 10 12 14 

L(fm) 



Figure 8. J/ip and ipf production cross sections versus L. The (parallel) lines correspond 
to exp(— paL), with pa = 0.088 fm~^. 
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ifjt cross sections. To cope with the range in the vertical scale, the branching ratios into 
muons were corrected for. The (S-U and Pb-Pb) L dependent a"^ values were obtained 
from t/'/DY using the dX\-E^ values for normalisation and assuming that DY does not 
depend on L (in other words, the integrated point is replaced by 5 points keeping the 
same average). The a"^ points were obtained from the ip' /if) values (with some art to 
combine the different centrality bins used in the ^/ip and ipi analyses |T^). 

We can see that the nuclear absorption model reproduces quite well all the i/ip points up 
to the most peripheral Pb-Pb collisions. The remaining Pb-Pb points reveal the existence 
of a new "suppression" mechanism, L not being any longer an appropriate parameter. 

The points reveal that after L ~ 6 fm the two charmonia states are already fully 
formed and start feeling in a different way the medium they cross. The nuclear matter 
path length alone is no longer enough to correctly parametrise the ipi behaviour, revealing 
that a dense hadronic medium is formed in S-U collisions. 

To understand whether the new Pb-Pb data on charmonia production is more than just 
a "hint" of QGP formation, accurate studies will have to be pursued, much beyond the 
"L framework". 

For instance, a correct understanding of the S-U data must consider that about 8 % 
of the observed ^/tp yield actually come from ipl decays. Since the ipi is significantly 
suppressed in S-U collisions, a certain fraction of the i/ip yield must also be suppressed. 

To conclude, table |^ collects useful information for further studies of i/ip production 
from p-A to S-U and Pb-Pb collisions, as a function of i?T- Notice, however, that the i?T 
scales are not the same in the NA38 and NA50 experiments, due to the different rapidity 
windows covered by the calorimeters: 1.7 < t] < 4.1 in NA38 and 1.1 < ?7 < 2.3 in NA50. 
The resolution of this non-trivial problem would require a few more pages. 

Table 7 

i?T dependence of J/ip yields: ip/DY and after normalisation to the integrated cross section 
and to the pp value. 







V^/DY 


S</(BA) 


</(BAO 




(GeV) 




(nb / nucleon) 


S-U 


All 


21.4±0.24 


l.OliO.lO 


0.48±0.07 




25 


25.2±0.69 


1.19±0.12 


0.56±0.09 




42 


22.8±0.56 


1.08±0.11 


0.51±0.08 




57 


21.0±0.48 


0.99±0.10 


0.47±0.07 




71 


20.2±0.44 


0.95±0.10 


0.45±0.07 




82 


19.2±0.44 


0.91±0.09 


0.43±0.07 


Pb-Pb 


All 


11.9±0.4 


0.59±0.04 


0.28±0.04 




35 


17.8±2.20 


0.88±0.13 


0.42±0.08 




59 


13.2±0.99 


0.65±0.07 


0.31±0.05 




88 


12.7±0.82 


0.63±0.06 


0.30±0.05 




120 


11.4±0.80 


0.57±0.06 


0.27±0.04 




149 


8.6±0.80 


0.43±0.05 


0.20±0.03 
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